Search for: All records

Fractures in Earth's critical zone influence groundwater flow and storage and promote chemical weathering. Fractured materials are difficult to characterize on large spatial scales because they contain fractures that span a range of sizes, have complex spatial distributions, and are often inaccessible. Therefore, geophysical characterizations of the critical zone depend on the scale of measurements and on the response of the medium to impulses at that scale. Using P-wave velocities collected at two scales, we show that seismic velocities in the fractured bedrock layer of the critical zone are scale-dependent. The smaller-scale velocities, derived from sonic logs with a dominant wavelength of ~0.3 m, show substantial vertical and lateral heterogeneity in the fractured rock, with sonic velocities varying by 2,000 m/s over short lateral distances (~20 m), indicating strong spatial variations in fracture density. In contrast, the larger-scale velocities, derived from seismic refraction surveys with a dominant wavelength of ~50 m, are notably slower than the sonic velocities (a difference of ~3,000 m/s) and lack lateral heterogeneity. We show that this discrepancy is a consequence of contrasting measurement scales between the two methods; in other words, the contrast is not an artifact but rather information—the signature of a fractured medium (weathered/fractured bedrock) when probed at vastly different scales. We explore the sample volumes of each measurement and show that surface refraction velocities provide reliable estimates of critical zone thickness but are relatively insensitive to lateral changes in fracture density at scales of a few tens of meters. At depth, converging refraction and sonic velocities likely indicate the top of unweathered bedrock, indicative of material with similar fracture density across scales. 

Abstract Erosion at Earth’s surface exposes underlying bedrock to climate-driven chemical and physical weathering, transforming it into a porous, ecosystem-sustaining substrate consisting of weathered bedrock, saprolite, and soil. Weathering in saprolite is typically quantified from bulk geochemistry assuming physical strain is negligible. However, modeling and measurements suggest that strain in saprolite may be common, and therefore anisovolumetric weathering may be widespread. To explore this possibility, we quantified the fraction of porosity produced by physical weathering, FPP, at three sites with differing climates in granitic bedrock of the Sierra Nevada, California, USA. We found that strain produces more porosity than chemical mass loss at each site, indicative of strongly anisovolumetric weathering. To expand the scope of our study, we quantified FPP using available volumetric strain and mass loss data from granitic sites spanning a broader range of climates and erosion rates. FPP in each case is ≥0.12, indicative of widespread anisovolumetric weathering. Multiple regression shows that differences in precipitation and erosion rate explain 94% of the variance in FPP and that >98% of Earth’s land surface has conditions that promote anisovolumetric weathering in granitic saprolite. Our work indicates that anisovolumetric weathering is the norm, rather than the exception, and highlights the importance of climate and erosion as drivers of subsurface physical weathering. 

In weathered bedrock aquifers, groundwater is stored in pores and fractures that open as rocks are exhumed and minerals interact with meteoric fluids. Little is known about this storage because geochemical and geophysical observations are limited to pits, boreholes, or outcrops or to inferences based on indirect measurements between these sites. We trained a rock physics model to borehole observations in a well-constrained ridge and valley landscape and then interpreted spatial variations in seismic refraction velocities. We discovered that P-wave velocities track where a porosity-generating reaction initiates in shale in three boreholes across the landscape. Specifically, velocities of 2.7 ± 0.2 km/s correspond with growth of porosity from dissolution of chlorite, the most reactive of the abundant minerals in the shale. In addition, sonic velocities are consistent with the presence of gas bubbles beneath the water table under valley and ridge. We attribute this gas largely to CO₂produced by 1) microbial respiration in soils as meteoric waters recharge into the subsurface and 2) the coupled carbonate dissolution and pyrite oxidation at depth in the rock. Bubbles may nucleate below the water table because waters depressurize as they flow from ridge to valley and because pores have dilated as the deep rock has been exhumed by erosion. Many of these observations are likely to also describe the weathering and flow path patterns in other headwater landscapes. Such combined geophysical and geochemical observations will help constrain models predicting flow, storage, and reaction of groundwater in bedrock systems.

 

As bedrock weathers to regolith – defined here as weathered rock, saprolite, and soil – porosity grows, guides fluid flow, and liberates nutrients from minerals. Though vital to terrestrial life, the processes that transform bedrock into soil are poorly understood, especially in deep regolith, where direct observations are difficult. A 65-m-deep borehole in the Calhoun Critical Zone Observatory, South Carolina, provides unusual access to a complete weathering profile in an Appalachian granitoid. Co-located geophysical and geochemical datasets in the borehole show a remarkably consistent picture of linked chemical and physical weathering processes, acting over a 38-m-thick regolith divided into three layers: soil; porous, highly weathered saprolite; and weathered, fractured bedrock. The data document that major minerals (plagioclase and biotite) commence to weather at 38 m depth, 20 m below the base of saprolite, in deep, weathered rock where physical, chemical and optical properties abruptly change. The transition from saprolite to weathered bedrock is more gradational, over a depth range of 11–18 m. Chemical weathering increases steadily upward in the weathered bedrock, with intervals of more intense weathering along fractures, documenting the combined influence of time, reactive fluid transport, and the opening of fractures as rock is exhumed and transformed near Earth’s surface.

 

We used seismic refraction to image the P‐wave velocity structure of a shale watershed experiencing regional compression in the Valley and Ridge Province (USA). From estimates showing strong compressional stress, we expected the depth to unweathered bedrock to mirror the hill‐valley‐hill topography (“bowtie pattern”) by analogy to seismic velocity patterns in crystalline bedrock in the North American Piedmont that also experience compression. Previous researchers used failure potentials calculated for strong compression in the Piedmont to suggest fractures are open deeper under hills than valleys to explain the “bowtie” pattern. Seismic images of the shale watershed, however, show little evidence of such a “bowtie.” Instead, they are consistent with weak (not strong) compression. This contradiction could be explained by the greater importance of infiltration‐driven weathering than fracturing in determining seismic velocities in shale compared to crystalline bedrock, or to local perturbations of the regional stress field due to lithology or structures.

 

The porous near‐surface layer of the Earth's crust – the critical zone – constitutes a vital reservoir of water for ecosystems, provides baseflow to streams, guides recharge to deep aquifers, filters contaminants from groundwater, and regulates the long‐term evolution of landscapes. Recent work suggests that the controls on regolith thickness include climate, tectonics, lithology, and vegetation. However, the relative paucity of observations of regolith structure and properties at landscape scales means that theoretical models of critical zone structure are incompletely tested. Here we present seismic refraction and electrical resistivity surveys that thoroughly characterize subsurface structure in a small catchment in the Santa Catalina Mountains, Arizona, USA, where slope‐aspect effects on regolith structure are expected based on differences in vegetation. Our results show a stark contrast in physical properties and inferred regolith thickness on opposing slopes, but in the opposite sense of that expected from environmental models and observed vegetation patterns. Although vegetation (as expressed by normalized difference vegetation index [NDVI]) is denser on the north‐facing slope, regolith on the south‐facing slope is four times thicker (as indicated by lower seismic velocities and resistivities). This contrast cannot be explained by variations in topographic stress or conventional hillslope morphology models. Instead, regolith thickness appears to be controlled by metamorphic foliation: regolith is thicker where foliation dips into the topography, and thinner where foliation is nearly parallel to the surface. We hypothesize that, in this catchment, hydraulic conductivity and infiltration capacity control weathering: infiltration is hindered and regolith is thin where foliation is parallel to the surface topography, whereas water infiltrates deeper and regolith is thicker where foliation intersects topography at a substantial angle. These results suggest that bedrock foliation, and perhaps by extension sedimentary layering, can control regolith thickness and must be accounted for in models of critical zone development. © 2020 John Wiley & Sons, Ltd.

 

The complex ecohydrological processes of rangelands can be studied through the framework of ecological sites (ESs) or hillslope‐scale soil–vegetation complexes. High‐quality hydrologic field investigations are needed to quantitatively link ES characteristics to hydrologic function. Geophysical tools are useful in this context because they provide valuable information about the subsurface at appropriate spatial scales. We conducted 20 field experiments in which we deployed time‐lapse electrical resistivity tomography (ERT), variable intensity rainfall simulation, ground‐penetrating radar (GPR), and seismic refraction, on hillslope plots at five different ESs within the Upper Crow Creek Watershed in south‐east Wyoming. Surface runoff was measured using a precalibrated flume. Infiltration data from the rainfall simulations, coupled with site‐specific resistivity–water content relationships and ERT datasets, were used to spatially and temporally track the progression of the wetting front. First‐order constraints on subsurface structure were made at each ES using the geophysical methods. Sites ranged from infiltrating 100% of applied rainfall to infiltrating less than 60%. Analysis of covariance results indicated significant differences in the rate of wetting front progression, ranging from 0.346 m min^−1/2for sites with a subsurface dominated by saprolitic material to 0.156 m min^−1/2for sites with a well‐developed soil profile. There was broad agreement in subsurface structure between the geophysical methods with GPR typically providing the most detail. Joint interpretation of the geophysics showed that subsurface features such as soil layer thickness and the location of subsurface obstructions such as granite corestones and material boundaries had a large effect on the rate of infiltration and subsurface flow processes. These features identified through the geophysics varied significantly by ES. By linking surface hydrologic information from the rainfall simulations with subsurface information provided by the geophysics, we can characterize the ES‐specific hydrologic response. Both surface and subsurface flow processes differed among sites and are directly linked to measured characteristics.